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Proceedings of former DERIVE & CAS-TI and ACDCA-Conferences 

You are invited to browse and to download lectures and workshops (all in df-format) of 
(almost) all earlier Conferences starting with Krems 1992 and ending with Malaga 2010: 

http ://rfdz.ph-noe . ac . at/ acdca/konferenzen.html 

Krems 1992, Krems 1993, Derive Days Diisseldorf 1995, Bonn 1996, Gettysburg 1998, 
Portoroz 2000, Liverpool 2000, Vienna 2002, Montreal 2004, Dresden 2006, Buffels- 
poort (South Africa) 2008, Malaga 2010 

Some proceedings have been published originally in printed form only. I scanned them 
all and interlinked them with the respective contents page. Many thanks to Peter Hof- 
bauer who set all files on the web. 

Happy Browsing and Happy Remembering for many of you - 
- and don’t forget: 

Attend TIME2012 in Tartu, Estonia 
July 11-14, 2012 
http://time2012.ut.ee/ 


Browse a library of over 2,400 educational videos at 

http://www.khanacademy.org 


Nur in deutscher Sprache (in German only): 

Franz Embachers Mathe-Online wurde erweitert und verbessert: 

http://www.mathe-online.at/neu.html 


DUG-Mitglied Johannes Zerbs ist sehr aktiv. Geme veroffentliche ich den Link zu seinen 
Publikationen und wiinsche weiterhin gutes Gelingen und viel Erfolg. 

MATHE CHILLED -ONLINE - ein Buch, eine Sammlung, ein "Ratselblock" von 
Ubungschilis(-kapiteln) nach den neuen Standards und Kompetenzen der Unterstufe 

http://www.iohanneszerbs.at/visionen/index.php?option=com docman&task=cat view& 
gid=4 1 &ltemid= 1 2 

Alle Infos mit der Buchbeschreibung unter Klick auf: 

— » www.johanneszerbs.at — » linker Meniipunkt „Online-Mathe-Ubungscenter“ — > 

—> UE-klasse-center — » UE-klasse-Worddokus 


Von der Homepage gelangt man auch zu einer mathematischen „Millionenshow“ und zu 
einem Mathe-Quiz, Josef 
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Liebe bUG-Mitglieder, 
eigentlich wollte ich in diesem Newsletter 
vor allem die Beitrage von bietmar Oertel, 
Piotr Trebisz und buncan Mcbougall wei- 
ter bzw. zu Ende fiihren. Ich war auch 
schon ziemlich weit mit den Ubersetzun- 
gen der deutschen Beitrage, aber dann 
erreichten mich einige Anfragen und 
mails, die schlieGlich liber weitere Kommu- 
nikation eine so heftige Eigendynamik 
entwickelten, dass ich sie nicht longer 
liegen lassen konnte. Nur Robert Setifs 
Anfrage musste auf den nachsten bNL 
verschoben werden. 

Ich mb'chte Sie besonders darauf hinwei- 
sen, dass nun die Tagungsbtinde aller bE- 
RIVE- und ACbCA-Konferenzen im pdf- 
Format verfugbar sind. Beachten Sie bit- 
te den Hinweis auf der Info-Seite. 

Ich freue mich, wieder einen TT-Nspire- 
Beitrag vorstellen zu konnen. Mit Version 
3 kann man sehr ordentlich arbeiten. Nun 
sind wirklich schon einige von uns lange 
gewiinschten bERIVE-Moglichkeiten imp- 
lementiert worden. 

Besonders spannend fur mich war die Be- 
schtiftigung mit der konformen Abbildung 
bzw. deren Umkehrung, zu der mich ein 
Brief von Klaus Ko'rner gebracht hat. 

Ich habe auch sehr wenig Ahnung iiber 
Tensore gehabt. bank Tom Fowler weiB ich 
nun ein wenig mehr daruber. Roger Folsom 
schreibt einiges liber Comparative Sta- 
tics" - wofiir ich keine passende deutsche 
Ubersetzung find en konnte. Ich habe aber 
ein paar sehr schone Unterlagen im Web 
gefunden, die zu einer weniger bekannten 
Anwendung der bifferentialrechung im 
wirtschaftlichen Bereich fiihren. 

Walter Klinger ist wieder zum Einsatz von 
bERIVE im Unterricht zuriickgekehrt und 
erlebt mit der Version 6 immer wieder 
Uberraschungen, die er gerne mit mir 
bespricht. Walter, vielen bank fur die 
anregenden Gesprache. 


bear bUG Members, 

Actually I intended to continue the con- 
tributions of bietmar Oertel, Piotr Tre- 
bisz and buncan Mcbougall in this news- 
letter. I was very busy translating the 
German papers but then I received some 
requests and emails which together with 
the resulting communication got a strong 
self dynamics and I could not leave them 
on my to-do-list any longer. Only Robert 
Setifs request must be left for the next 
issue. 

I'd like to inform you that the proceed- 
ings of all earlier bERIVE and ACbCA 
Conferences are available in pdf-format 
and can be downloaded. More info on the 
information page. 

I am very pleased about being able to 
present an TT-Nspire related article 
again. Version 3 allows comfortable and 
powerful working. Finally some of long 
desired bERIVE features have been im- 
plemented 

Extra exciting for me was dealing with 
the conformal mapping and its inverse 
provoked by Klaus Korner’s mail. 

I must admit that I knew nothing whatso- 
ever about tensors. Thanks Tom Fowler I 
now have at least a weak idea about these 
mathematical objects. Roger Folsom 
writes about "Comparative Statics" - I 
couldn't find an appropriate German 
translation. What I could find are some 
fine resources on the web. They are in- 
forming about a not so well known applica- 
tion of calculus in economics. 

Walter Klinger returned to use bERIVE in 
secondary level math teaching. He is com- 
ing across some surprises with bERIVE 6 
which he likes to discuss with me. Many 
thanks for the inspiring talks. 

With my best regards 

Mit herzlichen GriiBen ( Ai^Ai 


Download all ZWL-Derive- and TI-f!les from 

http : / / www . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 7/ User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


I II 

Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 


Next issue: December 2011 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encription, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 
Truth Tables on the TI, M. R. Phillips, USA 
Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 
Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ & Rob Gough, UK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier), CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Statistics with Tl-Nspire, G. Herweyers, BEL 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Rational Hooks, J. Lechner, AUT 
Cubus Simus, H. Ludwig, GER 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 
and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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On Quadratics and Conic Sections 1 

Wolfgang Propper, Nurnberg, Germany 
w.proepper@franken-online.de 


Nowadays conic sections have nearly disappeared from textbooks. Nevertheless it is an in- 
teresting field containing a lot of nontrivial algebra. 

1. (General 2-Variable) Quadratic Equations 

The intersection of a circular cone with the equation z 2 = x 2 + y 2 in R 3 with a plane 

ax + Py + yz = s yields a so called quadratic which only has terms in x 2 , y 2 , xy, x, y and 

a constant term. It looks like 


(1) ax 2 + 2b xy + cy 2 + 2dx + 2ey + f = 0. 

Due to its (geometrical) genesis its graph is called conic section or just conic. 

The left side of this equation is elegantly expressed by W J M W with M a symmetric ma- 


trix M = 


a 

b 

d ) 



b 

c 

e 

and W a vector W = 

y 

l d 

e 

n 


Iv 


. This construction leads to the factors 2 in 


the coefficients of the mixed and the linear terms [CS, 2.2], 


© the declarations: 


[l. 

Done 

term(xy):= 

((w(x^y )) T • ot • ))[ i l] Done 

eqn{xy)\=term[xy)=Q Done 

© for the examples we must declare a specific matrix M 

m: = 

a b d 
bee' 
d e f 

expand(e^/?(x 4 y)):© 2nd degree polynomial in- xand^ a x +2 b xy± 2 d x+c y +2 e y+f~ 0 

m:= 

a 0 d 
0 0 e '■ 
d e f 

term[xy)\Q a general expression for a parabola > a x +2 d x+2 e y+f 

| © now we become specific: 

m: = 

0 25 0 : PropFrac| ? 2 ^J+1:© <7/t ~ + V _1 

0 0 -225J 

m:= 

4 0 0 
0-9 0 
0 0 0 

(2x-3'_y)-(2x+3_y)=0 

.iactov\eqn\x,y))\G) a pair of straight line > 

m.= 

-1 0 C 

0 0 \ 

0 I 
2 

] y=x 2 +l 

^ :eqn(x,y}+x 2 +l :€> a parabola 

1 




11/99 


1 With this contribution is connected an extended Tl-Nspire document OnConicSections_NN.tns 

(NN is 20 or 30 depending on which release of Tl-Nspire CAS it should run. 20 is for release 2.x and 

30 for the latest release 3.x). When referring to this document in this contribution we denote this by 
[CS, Pr.Pg] where Pr means problem and Pg the page in this problem. 
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2. Classification of conic sections 

At first we will investigate the case b = d = e = 0. Equation (1 ) is reduced to 
(2) ax 2 + cy 2 + f = 0. 

In case a c f* 0 the classification is done easily by analyzing the signs of these coeffi- 
cients. 

If a, c and f have equal signs then the solution set is empty. Otherwise we have ellipses or 
hyperbolas, so called central curves with a centre point at the origin of the coordinate 
system and their half-axes on the coordinate axes. 


a 

c 

f 

Type 

± 

± 

+ 

Ellipse 

± 

+ 

± 

Hyperbola, symm. wrt y-axis 

+ 

± 

± 

Hyperbola, symm. wrt x-axis 

± 

± 

+ 

Empty set 


Table 1 


If however acf = 0, - i.e. at least one of the coefficients disappears -, then we have 3 in- 
stances: 

1 . Pair of straight lines (exactly one coefficient is zero): 

If f = 0 and a and c have different signs, the left side of equation (2) decomposes into 
two linear factors. This results in two symmetric straight lines through the origin. 

If a = 0 or c = 0 we get a pair of parallels with regard to the x- or the y-axis, if the signs 
of a or c and fare different. 

2. Double line (exactly two coefficients are zero): 

In case a and for c and f disappear then the x-axis or the y-axis is the only one solu- 
tion of the corresponding equation. 

3. Point: 

If f = 0 and a and c have equal signs the equation from above is only fulfilled by the 
origin. 

In all other cases there will be no more genuine intersection figure, i.e. the solution set is 
empty. Except if a = c = f = 0. Then equation (2) is fulfilled by the whole xy-plane. 


The following table shows this in a clear form: 
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a 

c 

f 

Type 

± 

+ 

0 

Pair of straight lines through the origin 

0 

± 

+ 

Pair of parallel lines with respect to x-axis 

± 

0 

+ 

Pair of parallel lines with respect to y-axis 

± 

0 

0 

Double line, the x-axis 

0 

± 

0 

Double line, the y-axis 

± 

± 

0 

Point, the origin 

± 

0 

± 

Empty set 

0 

± 

± 

Empty set 

0 

0 

± 

Empty set 

0 

0 

0 

Whole xy-plane 


Table 2 


The impact of coefficients b, d and e [CS, 2.3 and 2.4] 


If just b * 0 (and further on d = e = 0) then a rotation where x is replaced by 
xcos(#) - ysin(£) and y is replaced by xsin(6>) + ycos(#) can bring the coefficient of the 
mixed term to zero without introducing linear terms in x or y. This means that we have 
again one of the previous cases. 


System Matrix syst: 


a b d 
bee 
d e f_ 


© translation to remove linear terms ► 


koef^term[x+u,y+v )) 


Coefficient Matrix koef: 


System Matrix syst: 


2 2 

x xy y x y const 

a 2 b c 2-[au+bv+d) 2\b u+c v*-e ) au 2 +2u[bv+d)+cv 2 +2ev+f 

a b au+bvi-d 

b c bu+cv+e 

au+bv+d b u+c v+e u i/ 2 +2 i/ (i vtrf)+c v 2 +2'e vt/ 


Done 

solve 

IfcJSH 

be-cd , -[ae-bd] 

u= and v=— - 

ac-b 2 ac-b 2 


1 n, be-cd -( ae-bd ) 


tevtn\x+uy+v J|w— and v= 

, a c-b 2 a c-b 2 i 



Coefficient Matrix koef: 


System Matrix syst: 


2 2 
x xy y x y 


a 2b c 0 0 


ajef-e 2 )- 


b A f+2bde-cd z 


ac-b 


a b 
b c 


0 0 a\c f-e 2 )-b 2 ^2 b de-cd 2 


ac-b A 


© Result: A translation with u= ^ e c ^ and v= ^ ^ ° f lets the coefficients ofx and y disappear. 


ac-b * 


ac-b * 
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2/j 

The condition for removing the mixed term is tan(20) = . (In case a = c the conic sec- 

a-c 

tion is a circle which is invariant with respect to rotations around its centre.) 

Likewise (at an arbitrary b) and d, e * 0 the translation x := x + u and y\=y+v may re- 

b ' g — c ‘ d 

move the coefficients of the linear terms in x and y. The conditions are u = and 

y ac-b 2 

v = . The coefficients of x 2 , xy and y 2 remain unchanged. 

ac-b 


© Here we determine a rotation that the coefficient of x*y disappear* 


© first we assign syst (from page 3) to m: 


=syst 


a b 
b c 


0 0 a\f±e^± 


a-c-b 


© Statement to calculate the angle of rotation* 


koeJ^term[x- cos(0)- i y sin(0),x-sin(0)+y-cos(0))) 


:oefficient Matrix koef: 
2 


*-y y 

2 +(2 A cos(0)+c sin(0)) sin(d) -2-(a sin(0) cos(^)-i' (2 (cos(0)) 2 -l)-c sin(0) cos(^)) a- (sin(^)) 2 — (2 b sin(0)-cccs(0))cos{0) 0 0 a ^ c f e ^ b f +2 ^ a 


x y const 

.2U2., 


iystem Matrix syst: 


a(cos{0)) 2 +(2Z>cos(0)+csin(0))'sin(0) -<7sin(0)cos(0)+b’ (2(cos(0)) 2 -l)t-c- 

o sin(^) cos{^)+Z>- (2 (cos(0)) 2 -l)+c sin(0) cos(0) t7-(sin(^)) 2 -(2- isin(0)-ccos (0))cos(0) 


0 

0 

a{cf-e 2 )-b 2 f+2b de-c d 2 


ac-b 


© now we have to calculate 0 that the mixed term disappears* 


solve(sys/[ 1, 2]=0, d) 


'2 b] . 
l+nl-n 

iQ-c} 


© result: The condition forQ is: tan(20)=- 


2b 


© because m was not changed we can now take the statement from above and substitute (* 


^o^/erw(x cos(^)->' sin(0),x sin(0)+>' cos(0))|^=-^- tan -1 j 


There is a lot of algebra in the Calculator Window, isn’t it? 


As one easily can see, these conditions are only valid if a-c - b 2 ^ 0. This case has not 
been discussed up to now. Here we have a parabola, as we will see in the next section. 


The parabola case [CS, 2.5] 

If ac - b 2 = 0 and a^Owe can omit c in equation (1 ) and we get 

9 b 2 7 

(1 ') ax 2 +2bx- y + — y 2 + 2d • x + 2e • y + f = 0 

a 

This equation can be brought to the form 

(3) a'x 2 + 2e'y = 0 

with a' and e' only depending on the coefficients of equation (1 ’). (If in (T) the term with x 2 
should be omitted then (3) takes the form c'y 2 + 2 d'x = 0. The calculations are principally 
the same.) 
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In both cases we have parabolas, so called non-central curves. They are symmetric with 
respect to the y- or x-axis and they are passing the origin. 


b 

Applying a rotation by angle #with tan(#) = — we can remove the mixed term and the 

a 

term which is quadratic in y in equation (1') as well. So it takes the form 

a’x 2 + 2d'x + 2e'y + f = 0. 


d' d' 2 -a' f' 

And after a translation x := x + u and y := y + v with u = — and v = we finally 

a' 2a' e' 


yield equation (3). 



Conclusion 

We can constrain the investigation of equation (1) to two cases: 

It is sufficient that either b = d = e = Oor that b = c= d= f=0 (or a = b = e = f). The first 
case yields ellipses and hyperbolas and their degenerations (straight lines or point). In the 
other case the conic section is a parabola. 

3. Classification using a matrix M [CS, 3.1-3] 

By an operation H J M H, where matrix H is either a rotation matrix 


R = 

'cos (0) -sin(0) (f 
sin(0) cos(Q) 0 

or a translation matrix T = 

(1 0 u) 
0 1 v 

, matrix M = 

r a b d ' 

bee 


l 0 0 lj 


to 0 lj 


[d e f) 
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can be transformed to 

fa' 

0 

0 0^ 
c' 0 

, and, in case a c - b 2 = 0, to 

fa' 0 0^1 

0 0 e' 


1 ° 

0 f’J 


[0 e' 0) 


tively to 


0 0 d 

0 c' 0 
d' 0 0 


'A 


(or respec- 


). As the kind of a conic section is neither changed by a translation 


nor by a rotation, we can use matrices for the classifications provided in paragraph 2. 

For doing this we introduce new quantities: A = det (M), 5 = det ^ j and S = a + c. 
They are invariant with respect to the rotations and translations from above. So we get: 


8 

A* 0 

A = 0 

>0 

zl S<0: Ellipse 

Point 

Empty set 

<0 

Hyperbola 

Intersecting pair of straight lines 

=0 

Parabola 

Pair of parallel straight lines 


Table 3 



Again a lot of algebra! 


4. Representation of Conic Sections 

In this section some constructional ideas for conic sections are presented. All these ideas 
use the ‘zeros’ command from Tl-Nspire CAS. (So they can not be used in the nonCAS 
version of Tl-Nspire.) 
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If p(x,y) is a 2 nd degree polynomial in x and y then zeros(p(x,y), y) yields a list of (root-) 
expressions in x. The G&G application interprets this list as a family of functions in x with 
two elements maximum. So a definition f1(x):=zeros(p(x,y),y) in G&G plots a graph with 


two colors and labels f1_1 and f1_2. 


Many thanks to Philippe Fortin (France) for providing most of the ideas realized in the following 
paragraphs. 


4.1 The Conic Section Generator [CS, 4] 

The classification of conic sections from paragraph 2 and 3 can easily be studied with 
the ConicSectionGenerator (CSG). 


The workspace is divided into two G&G windows, a (smaller) Geometry and a (wider) 
Graphs window. In the left window we insert 6 sliders for variables a through f having 
a range from -10 to 10, which can easily be changed. The step size is 0.1 for a, b, c, 
and 0.5 for d, e, f and may also be changed as needed. 

The right window has the definition 

f1(x):=zeros(a x 2 + 2 bxy + cy 2 +2 dx + 2 ey +f= 0, y), (in function mode) which plots 
the (maximum) two components of the conic section. 


In addition the determinants A = det 


a 

b 

b 

c 

d) 

e 

, S= det(f 


[d 

e 

f) 

l c 

c ) 


and the sum S = a + c 


are shown in the upper left corner. 


In the upper right corner we see angle of inclination of the main axis and the coordi- 
nates of the centre point. 


By modifying parameters a through f ellipses, hyperbolas, parabolas and straight lines 
can be created and can be compared with the summaries given in tables 1, 2 and 3. 
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This is the screen after moving some sliders, i.e. changing some coefficients: 



Below are the screens from above as they are presented on the handheld Nspire: 




4.2 The Focus Feature [CS, 5] 

For ellipses the sum of the focus rays, i.e. the line segments connecting one point on 
the ellipse with the two focus points is constant. That means d(F^X) + d(F 2 X) = k 
(= const.). For a hyperbola this condition is \d(F 1 X) - d(F 2 X)\ = k. Both expressions on 
the left side are containing roots. One can eliminate these roots by quite elaborate al- 
gebraic transformations [see CS, 5.2] in order to obtain a rather complicated expres- 
sion in x and / which depends on the focus point coordinates fn,f ] 2 , f 21 , f 22 , and the 
constant k. 

The value of k can be varied by a slider. For the focus point coordinates we store the 
coordinates of two arbitrary points in fl 1 , fl 2 and f21 , f22 respectively. Now the graph 

can be plotted, using again the zeros function. By dragging the foci and changing /ewe 
can create ellipses and hyperbolas [see CS 5.3]. 

In general it would also be possible to insert the focus point condition directly in G&G, 
that means to define fl (x) = zeros(abs(norm([x;y]-[f11;f12])-norm([x;y]-[f21;f22]))-k,y). 
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If previously the coordinates and k have been created like mentioned above we should 
get the graph of a conic section. However, it takes quite a long time (depending on the 
performance of the computer) to receive a result. And if the foci are dragged around or 
k is changed Nspire is blocked for several minutes and the dynamics of the application 
gets lost. 


A final remark: As for this example there are not so many objects on the screen, it can 
also run on a Tl-Nspire CAS Handheld device (with reasonable visibility). The compli- 
cated calculations from [CS, 5.2] are omitted and just the result is inserted. The file 

Focus_NN.tns, where again NN=20 indicates a document for Nspire version 2.x and 
NN=30 is for version 3.x is among the files connected with this DNL. 



Ql 5.1 1 5.2 II 5.3 

► *OnConicSe..._30 ^ 

ik = 5. 

! "■ 

Art cfes -Hyper 


i i i i i ! 

2 °.j 

A,34.4,-2.9) 

jJ 



0 



2 



f m F2 

# 


1 

1 


Two handheld-screenshots of FOCUS 


4.3 The Five Point Feature [CS, 6] 

The quadratic a-x 2 + 2 bx-y + c-y 2 + 2 d-x + 2 e-y + f = 0 in x and y contains 6 coeffi- 
cients. With the coordinates of five different points inserted for x and y we get five 
equations. An equation like a + 2b + c = 1 avoids that a, b and c are zero at a time and 
we call it normalization condition. 

This system yields a unique solution for a through f m . 

For the implementation we define a function fivepoint(P1, P2, P3, P4, P5) where 

PI,..., P5 are 2-dimensional column vectors with components pi,..., p5 and ql q5. 

This function solves the system mentioned above and assigns the solution to 

a-x 2 + 2 b-x-y + c-y 2 + 2 d-x + 2 e-y + f. 

We only have to draw five different points in a G&G window, assign their coordinates 
to pi,..., p5 and ql,..., q5 and define f1(x):=zeros(fivepoint(P1,P2,P3,P4,P5),y). A 

conic section through those five points appears. Now we can grab one or another 
point and drag it around observing what is happening. 

I have provided a stand alone document to run on handhelds. It is called 
FivePoint_NN.tns and consists of only two pages: a short introduction and the con- 
struction in G&G. The function fivepoint is contained in the document but no longer di- 
rectly visible. 

Two screenshots of fivepoint are given on the next page. 
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[1] Additional Comments of the Editor: 

It might be a good question for students, how to explain the fact that there are 6 unknown coefficients 
in the quadratic but only five points are sufficient to define a conic? 

d 2 b 

Dividing (1 ) by f yields six 2 + b'xy + c'y 2 + d'x + e'y + 1 = 0 with a' = —,b' = ... . Now it is easy to see 

that 5 points are sufficient. (One could also divide by a, ...). 


Tl-Nspire has implemented a long desired DERIVE- feature: Background pictures. 

I load a picture of a bridge (jpg, bmp, ...) not far from my home and would like to analyze the form of 
the arc. Ask students (and teachers!!) and there first answer will be: This is a parabola. 

Now drag the five points as accurate as possible on one of the arcs and you will find a wonderful 
curve. Go back to the Calculator App and ask for fl(x). As Wolfgang mentioned above you will find two 
branches, take one of them, do a little algebra and get the quadratic. Then follow the recipe given in 
CS for classifying and analyzing the curve. Or you write a program which does the job for you! 

You could also use the CSG (page 10) to model the bridge (or other pictures). This is a little bit more 
sophisticated but the students might find out the impact of the various parameters. 

Finally, there is a German version of this contribution at 

http://www.ti-unterrichtsmaterialien.net/imqserv. php?id=Pr6pper.pdf 
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Let’s do a little algebra to find the quadratic in its implicit form: 


/**) 


-37.4498 (Jx 2 +10.4995 x+34.7117 -0.910278'(x+6.55427)), 37. 4498 (,/x 2 +10.4995 x+ 34.7117 +0.910278 (x+6.55427)) 


'=//W[i] 


'=-37.4498'(Jx 2 +10.4995x+34.7117-0.910278'(x+6.55427)J 


=-37. 449809063834' (Jx 2 + 10. 499493646781 x+34. 71 1708292528 -0. 91027781213277' (x+6. 5542685029883)) 

-37.4498 

-0. 026702 y=l . • (Jx 2 +10.4995 x+ 34.7117 -0. 910278- (x+6. 55427)) 


-0.026702412295927-^=1. 0000002420262- (Jx 2 +10. 499493646781 x+34. 71 1708292528 -0.91027781213277- (x+6. 5542685029883)))+0. 910278' (x+6. 55427) 

0. 910278 x-0.026702>+5. 96621=1. Jx 2 +10.4995 x+34.7117 -3. 24439E-8-X+0. 000001 


(o. 910278 x-0. 026702412295927 -_y+5. 96620778706=1 . 0000002420262' Jx 2 +10. 499493646781 x+34. 71 1708292528 ) 2 

0. 828606- (x-0. 029334- (y-223.433)) 2 =l . (x 2 +10. 4995 x+34. 71 17) 


expand(lefl(o. 828606037284- (x-0. 029334348732944- (y-223.43328838384)) 2 =l. 0000004840525' (x 2 +10. 499493646781 x+34. 71 1708292528))-right(o.82860603 > 


-0. 171394'X 2 -0.048613-X i y+0.362317-x+0.000713- < y 2 -0.318624 , - y+0. 88391=0 


The next contribution will do this job for you! 
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A Conics Explorer and a Conics Trainer 

Josef Bohm, Wiirmla, Austria 


As promised on page 13 here is the analysis of the “bridge - conic”. 

2 2 

task := - 0.171394-x - 0.048613-x.y + 0.362317-x + 0.000713-y - 0.318624-y + 0.88391 = 0 


ans 


Type 

Center 

Real axis 2a 

Imaginary axis 2b 

Real Vertices 





-10.99451200 85.94367265 

Hyperbola 

[-5.249788684, 44.47126414] 

83.73679034 

12.79439536 

0.4949346356 2.998855624 


Asymptotes 

0. 09253588869- t - 5.249788684 6 . 620111512 • t + 44.47126414 

1. 888548332- t - 5.249788684 44.47126414 - 6 . 345693384 . t 


It seems to be a hyperbola. I say “it seems to be” because unlike to translation and rotation the conics 
can change their type after applying a perspective projection which is given by taking a picture of the 
object with a camera. So it could be indeed a parabola! ! 

In fact, the program would not give the answer parabola - even if it would be one - because of 
numerical reasons. As Wolfgang explained on page 6 the discriminant must become zero - and I am 
quite sure that the pixel coordinates would not allow this. 

This example was the introduction. 

The Conics Explorer is the core part of my Conics Trainer: 

Load conics . mth as a Utility file and simply Simplify start 

#1: L0AD(D: \DOKU5\DNLs\DNLll\MTHS3\coni cs . mth) 

#2: start 

Conics Trainer - Trainer fur Kegel schnitte 

coni = Center in origin un rotated - Kegel schnitt in Ursprungslage 
con2 = shifted origin, unrotated - verschobener Kegelschnitt 
con 3 = conic shifted and rotated - Kegelschnitt in al 1 genie iner Lage 
con = random choice from above - Zufal lsaufgabe aus den obi gen 

You can enter your own conic as task:= ... 


ans = gives the analysis; ant = liefert die Analyse 
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Let’s assume that your students should train the conics of type 2 (shifted origin and not rotated): 
con 2= delivers a randomly generated quadratic. 

The students do their investigation with or without technology and then they can check their solution. 


2 

#4: con2 = (9-x - 16-y - 7-y = -30) 




Type 

Vertex 


Parameter 

Focus 


#5: 

ans = 

Parabola 

1969 

7 

9 

59 

7 



576 

32 . 

32 

18 

32 . 



Finally they can plot the curve together with all the important details. 



con 3= gives the most complicated cases - shifted and rotated. 

2 2 

#6: con3 = (12-x + 9-x-y + S-x + 12-y - 17-y = -16) 

#7 : ans 


#8: 


#9: 

#10: 


Typ 

Mittelpkt 

Hauptachse 

Nebenachse 

Scheitel 

Type 

Center 

Major axis 

Minor axis 

Vertices 

. Ellipse 

[-0 . 6969696969 , 0 . 9696969696] 

imaginary 

imaginary 

imaginary vertices . 

con3 = (8- 

2 2 
x - 5-x-y - 13-x + 15-y + y = 

21) 






Typ Mittelpkt Hauptachse Nebenachse 

Type Center Major axis Minor axis 


Scheitel 

Vertices 


# 11 : 


Ellipse [0.8461538461, 0.1076923076] 3.833364128 2.587418828 


2 . 6714015 67 0 . 6926204161 

-0 . 9790938750 -0 . 4772358007 
0 . 45134293 06 1.339685 972 


L 1.240964761 -1.124301356 J J 
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The quadratic in expression #6 turns out as an imaginary conic, but in expression #9 we find a rotated 
and shifted ellipse. 



The program code: 

start_ (dummy) := 

Prog 

dummy := RANDOM (0) 

[x :e Complex, y :e Complex] 

DISPLAYC' 1 Conics Trainer - Trainer fur Kegel schni tte ri ) 

DISPLAY( M conl = Center in origin unrotated - Kegelschmtt in Ursprungslage' 1 ) 
#1: DISPLAY( M con2 = shifted origin, unrotated - verschobener Kegelschmtt' 1 ) 

DISPLAYC M con3 = conic shifted and rotated - Kegelschmtt in allgerneiner Lage M ) 
DISPLAY( M con = random choice from above - Zufallsaufgabe aus den obigen M ) 
DISPLAYC' 1 M ) 

DISPLAYC M ans = gives the analysis; ant = liefert die Analyse' 1 ) 

DISPLAYC" M ) 

#2: start := start_C) 


#3: task := 

coni cl (a., b, c, rd) := 

Prog 

rd := RANDOM (8) 

Loop 

a := 20 - RANDOM (41) 
b := 20 - RANDOM (41) 

#4: f := 30 - RANDOM (61) 

If f-a-b t 0 
Prog 

task := [a - x A 2 + b-y A 2 - f = 0, a-x A 2 + b-y = 0, a-y A 2 + b-x = 0] 
task := IF(rd < 6, taskj.1, IF(rd < 7, task|2, task; 3)) 

RETURN task 
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conic2(a, b, c, d, e, rd) := 

Prog 

rd :: RANDOM (8) 
a :: 20 - RANDOM (41) 
b := 20 - RANDOM (41) 

#5: c := 20 - RANDOM (41) 

d := 20 - RANDOM (41) 
e :: 30 - RANDOM (61) 

task := [a*x A 2 + b-y A 2 + c-x + d-y + e = 0, a-x A 2 + c-x + d-y + e = 0, b-y A 2 + c*x + d-y + e = 0] 
task := IF(rd < 6, taskjl, IF(rd < 7, task|2, task|3)) 

RETURN task 

comc3(a, b, c, d, e, f) := 

Prog 

a := 20 - RANDOM (41) 
b := 20 - RANDOM (41) 
c := 20 - RANDOM (41) 

#6: d := 20 - RANDOM (41) 

e :: 20 - RANDOM (61) 
f := 30 - RANDOM (61) 

task := a-x A 2 + b-y A 2 + c-x + d-y + e-x-y + f = 0 
RETURN EXP AND( task) 


#7: 

coni 

:= coni cl () 

#8: 

con2 

:= coni' c2() 

#9: 

con3 

:= conic3() 

#10: 

con : 

= [coni, con2, con3] 

RANDOM (3) + 1 


Expression #1 1 (not printed here) is the main function ch_coni c() . It investigates the random 
conic con from #10 which is also stored as global variable task and stores the results in ans and ant 
(for the German users!). 

#12: ans := ch_coni c(task) 

#13: ant := ch_coni'c(task) 


Let’s try finding a hyperbola: 

2 2 

con3 = (18 -x +34-x-y-13*x+6-y +17*y= -12) 

ans 

Type Center Real axis 2a Imaginary axis 2b 

Hyperbola [-1.013812154, 1.455801104] 4.532953869 2.030942852 


Real Vertices 

-2 . 322866802 3 . 306015585 

0 . 2952424933 -0 . 3944133754 


+ 1.455801104 


Asymptotes 


- 0. 4727708474* t - 1.013812154 
2 . 105449700 * t - 1.013812154 


2 . 399597608 * t 
1.455801104 - 



1. 244452697- 1 
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The next task is not so easy. 

Model the arc as a conic using the slider bars: 



You can find conics everywhere. I used the sliders to model the first three ellipses. Three more are 
remaining waiting to be plotted! 
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Textbooks and the arctanh 

An email-discussion 


From: Michel Beaudin 

To: Albert Rich, David Jeffrey, David Stoutemyer, Josef Boehm 
Subject: Textbooks and the arctanh 

Hello Albert and my friends - it is always exciting for me when I am writing an email to the 
"fathers of Derive", to my best ROC friend ("ROC" stands for Rest of Canada: Quebecers love 
this word) and to my best friend in Europe. 

My question is: do you know why textbooks (at least calculus ones) are so quiet about the 
arctanh function? I understand that arctanh(x) has the open interval -1 < x < 1 as domain 
and this can cause a problem for students when they will be integrating 1/(1 - x A 2) if the 
bounds of their definite integral are, for example, 3 to 5. But it is quite easier for them to 
use the antiderivative arctanh(x) instead of 2 logs (or one log with quotients) because of the 
symmetry with the antiderivative arctan(x) for 1/(1 + x A 2). Moreover, when textbooks will 
start to use arctanh(x) for the antiderivative of l/(l-x A 2), it will be natural to perform the 
integral of sec(x): sec(x) = l/cos(x) = cos(x)/cos A 2(x) = cos(x)/(l - sin A 2(x)), so the sub- 
stitution u = sin(x) leads to l/(l-u A 2) and that will be a good reason to use arctanh(sin(x)) 
instead of log (abs(sec(x) + tan(x))) - what textbooks are still using - as an antiderivative 
of sec(x). 

Looking back to an old copy of the DNL journal, David Stoutemyer wrote: "Calculus students 
are introduced to complex numbers in high-school algebra, and they are taught that 
antiderivatives can differ by arbitrary constants. Why not exploit logarithms of negative 
numbers to reinforce both ideas? This is not to say that we should generalize all elementary 
calculus to complex numbers — only that if a complex number is most correct and natural in 
a few places, exploit it rather than fight it." (DNL #26, page 5). 

Things do not change as fast as I would like! 

Regards, 

Michel 


Von: Albert Rich 

An: Michel.Beaudin, David Jeffrey, David Stoutemyer, Josef Boehm 
Betreff: Re: Textbooks and the arctanh 

Hello Michel, 

Not being an educator, I am probably least qualified of your addressees to know “why textbooks are 
so quiet about the arctanh function?” Or, more generally, why the inverse trigonometric functions get 
so much more "ink” than the inverse hyperbolic ones? This is despite the undeniable symmetry be- 
tween these two sets of functions. This symmetry is evidenced by the identical, parallel treatments 
they receive in Chapter 4 of the recently released NIST Handbook of Mathematical Functions, which 
is the authoritative successor to Abramowitz and Stegun. 
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So good luck on your quest to make the inverse hyperbolic functions (and hyperbolic functions?) first 
class objects of study. 

The father of Rubi, 

Albert 


De : David Jeffrey 

Objet : Re: Textbooks and the arctanh 

My guess is that trigonometric functions are introduced early in high school because of the 
geometric application. 

Hyperbolic functions, on the other hand, have to wait until the exponential function has been 
covered, and thereafter they never catch up to trig functions in the mathematical mind. 

One advantage of arctanh(sin x) - or invtanh(sin x) as I prefer to write it - is that it is always a 
real value, unlike ln(sec x + tan x). 

For the same reason arctanh( (x A 2-1)/(x A 2 + 1) ) is a better integral of 1/x than ln(x) but I don't 
think it will catch on. 

DJJ 


From: Michel Beaudin 

Subject: RE: Textbooks and the arctanh 

David, if the first goal when choosing an antiderivative is to have a real value, then I "prefer" 
ln(abs(x)) instead of arctanh((x A 2-l)/(x A 2+l)) for the integral of 1/x (these 2 expressions 
are identical). In the same idea, ln(abs(sec x + tan x)) is always real, so those who are writ- 
ings calculus textbooks will continue to use it instead of arctanh(sin x) for the integral of 
sec x. 

Michel 


Von: David R Stoutemyer 

Betreff: Re: RE: Textbooks and the arctanh 

Hi Michel, 

I can remain quiet no longer. 

The absolute value: 

1. unnecessarily makes the antiderivative incorrect for non-real endpoints (and some people 
will eventually use non-real endpoints not realizing the incorrectness), 
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2. completely hides the non-integrable singularity if someone integrates from a negative to a 
positive value. (With ln(x) you get a some warning from the i pi term.) 

3. makes it very difficult to simplify subsequent expressions that use the antiderivative -- for 
examples, mixtures of ln(x) and ln(abs(x)). As another example, iterated integrals become 
very difficult if an inner integral generates a ln(abs(x)). 

With all of its extra closure, the complex domain is so much easier than the real domain. I 
don't know why we put these extra hurdle constraints on students who have been exposed to 
complex numbers in earlier courses. 

-- aloha, david s. 


Von: Michel. Beaudin 

Betreff:* RE: RE: Textbooks and the arctanh 
Hi David, 

You are absolutely right and I will add more! 

After their calculus course, my students will take the ODE course. They will learn how to 
solve by hand "easy" ODEs and also they will learn how to use the "desolve" command of 
their CAS calculator (Voyage 200 and, starting in September, Nspire CX CAS). Those who 
have set the "complex format" to "Real" will sometimes get a "Non-real result" message on 
their screen when trying to solve some ODE. With complex format set to "Rectangular", we 
get a solution involving the number "i" but, choosing the appropriate domain - and this is a 
nice opportunity to recall the existence and uniqueness of solution ~, the "i" will "disappear" 
and students will obtain a real solution! 

In my talk at ACA 2007 ("Using ln(abs(x)) as an Antiderivative for 1/x is a Bad Choice!"), I 
used the separable ODE y ' = 4*y/(x A 2-9), y(0) = 5. In order to remove the imaginary unit 
"i" from the solution found by the device (and also Derive, Maple, ...), students have to use 
the fact that we are looking for a solution defined in the strip - 3 < x <3 because the initial 
condition is at the point (0, 5) : according to the existence and uniqueness theorem, this 
ODE has a unique solution defined for -3 < x <3 because both f and dif(f, y ) are continuous 
in some rectangle centred at the point (0, 5) with -3 < x < 3 (here f stands for 4*y/(x A 2 - 
9)). And restricting the domain gives this unique and real solution! 

Consequently, we are doing more and not less mathematics if we avoid absolute values in 
the integral of 1/x. 

Best regards, 

Michel 
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On 11-05-17 03:00 PM, Josef Bohm wrote: 

Dear friends, 

I am very fascinated by your high level discussion. Unfortunately I can not contrib- 
ute because I never had to do with hyperbolic functions and their inverse. BTW, in 
German the inverse hyperbolic functions are called “Area- functions”, eg. Area sinus 
hyperbolicus, etc. I never taught on a technical school where the hyperbolic func- 
tions might be part of the curriculum. 

What I can contribute is the fact that in Austria secans and cosecans are widely 
unknown and unused. I am quite sure if I’d ask some math teachers about secans 
and cosecans - they would have never heard about them. They don’t appear in our 
textbooks. I never heard about sec and cosec during my math study on the Techni- 
cal University ... 

Many thanks and I will follow your discussion with high interest - and I will not 
hesitate including it into the next DERIVE newsletter. It is a nice coincidence that I 
am just revising DNL#26 with David’s paper “In x or In | x | ?” 

Best regards to you all, 

Josef 


From: David Jeffrey 

Subject: Re: AW: RE: Textbooks and the arctanh 
Dear All, 

My Russian wife has been teaching North American calculus for the first time, and also com- 
mented that she was not taught secant and cosecant. I guess in Europe you are taught to write 
l+tan^ = 1/cos ^2 instead of our 1 +tan' N 2=sec /v 2. 

We always boast that mathematics is a universal language, but maybe not quite as universal as 
we hope. 

Some years ago, James Davenport was working on the open math project and reminded us of 
even more obscure functions such as versin(x)= 1 -cos(x) See 

http://en.wikipedia.0rg/wiki/File:Circle-triq6.svq 


for an attractive diagram. 
David J 
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Von: David R Stoutemyer 

Betreff: Re: AW: RE: Textbooks and the arctanh 

Tables of esc, sec and cot usefully replaced a division with a multiplication in the days of 
manual computation. However, those days are long gone, and most students quickly forget 
even whether esc is the reciprocal of cos or sin. Therefore we should stop wasting valuable 
curriculum time torturing students with these secondary trigonometric functions and their hy- 
perbolic counterparts, just as we no longer torture students with other trivially related trigo- 
nometric functions such the haversine, which was convenient for manual celestial navigation 
computations. We don't have a special name for the reciprocal of ln(...), abs(...), etc. 

Also, I think that hyperbolic functions are way less important than trigonometric, because 
exponentials are often a reasonable alternative to hyperbolics, whereas complex exponen- 
tials are a less attractive alternative to sinusoids. (Yes, even I prefer results absent i when i 
offers no strong advantage.) 

-- aloha, david s. 


Von: Michel. Beaudin@etsmtl.ca 
Gesendet: Mittwoch, 18. Mai 2011 02:22 
An: nojo.boehm@pgv.at 
Betreff: RE : RE: Textbooks and the arctanh 

Dear Josef, this email is sent only to you because, some months ago, I sent my "new integral 
table" to David Jeffrey and Albert Rich but not to you (David Stoutemyer also received it). In 
fact, this "high level discussion" started last November when I was teaching single variable 
calculus (see the abstract of my Derive session) and when I started to be fascinated by Rubi! 
So, here is the table of integrals that I would like to use and also here is the table actually 
used at ETS in the single variable calculus course. If you can't make it to Houston, this table 
of integrals will give you an idea of what I will be talking about in my lecture! 

Best regards, 

Michel 



Do you know this solid? 
Find more on page 38! 
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This is Michel Beaudin’s very individual table of integrals (which he is finding to be 
sufficient): 
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An email came in, which led to an interesting application (translation of the 
most important parts of the mail next page). 

Sehr geehrter Herr Bohm, 

durch Zufall bin ich auf Ihre e-mail Adresse gestoften und mochte Sie, wenn es fur Sie keine 
allzu grofte Miihe macht, um eine mathematische Hilfe bitten. Ich habe gehort, dass die Lo- 
sung meines Problems mit DERIVE moglich ware, doch kann ich selbst damit nicht umge- 
hen. 

Einige Satze zu meiner Person: Ich bin 1939 geboren, habe in Miinchen an derTU Maschi- 
nenbau studied, wardann mein ganzes Berufsleben in einer Turbinenfabrik in Niirnberg be- 
schaftigt. In relativ jungen Jahren wurde ich zum Leiter der technischen Berechnung und der 
technischen EDV-Abteilung gemacht und hatte mit Thermodynamik, Stromungsmechanik, 
Festigkeitslehre, Schwingungslehre, Materialfragen, Lagerauslegung, Wellenschwingungen 
und Wuchttechnik zu tun. Mein ganzes Berufsleben hat mir viel SpaR> gemacht. Als "techni- 
sches Gewissen" unserer Vertriebsleute hatte ich die Moglichkeit fast die ganze Welt zu be- 
reisen doch gesehen hatte ich dabei relativ wenig weil nie Zeit war. 

Nun im Ruhestand denke ich manchmal uber Probleme nach die mich wahrend meines Be- 
rufslebens irgendwie begleitet haben, doch deren Losung aus Zeitmangel immer verschoben 
wurde. 

Es geht um das Wissensgebiet "konforme Abbildungen". Wenn Sie sich auch damit ausken- 
nen, dann iiberlesen Sie bitte meine nachfolgenden kurzen einfiihrenden Worte. 

Mittels konformer Abbildung kann man einen Kreis in eine tragfliigelformige Figur verwan- 
deln. Diese Figur sieht dann auch so ahnlich aus wie ein Schaufelprofil. Meine Idee ist es 
nun ein gegebenes Schaufelprofil zuriick zu transformieren und zu sehen, ob die Rucktrans- 
formation eine "schone, ebenmafiige" Figur ergibt oder so etwas wie eine "Kartoffel". In die- 
sem Fall wurde ich daraus schlieften, dass dieses vorgegebene Profil stromungstechnisch 
nicht optimal ist. 

Gegeben sind die Koordinatenwerte eines Kreises mit vielen x =... und y = ... als Funktion 
von Radius und den Mittelpunktskoordinaten. 

Mit den Formeln u = x * (1 + k/(x 2 + y 2 )) und v = y * (1 - k/(x 2 + y 2 )) ergibt sich als transfor- 
mierte Figur eine Schaufelprofil in einem Koordinatensystem u,v. 

Mein Problem liegt nun darin, aus diesen beiden vorgenannten Formeln u = ... und v = ... die 
Werte x und y als Funktion von u und v darzustellen. Ich habe mir schon die Finger krumm 
gerechnet und bin zu keiner Losung gekommen. 

Man sagte mir allerdings, dass eine solche Formelumkehr mit DERIVE moglich ware. 

Da Sie sich offensichtlich mit diesem System auskennen, wage ich es, Sie hier um eine fach- 
liche Hilfe zu bitten. 

Ich wurde mich sehr tiber eine Antwort freuen, gegebenenfalls auch auf einen Hinweis wie 
ich hier weiterkommen konnte. 

Viele Griifte aus Nurnberg 

Klaus Korner 
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Klaus Korner’s problem - Inverse Conformal mappings 

... My problem deals with the field of knowledge “conformal mappings 

Applying a conformal mapping it is possible to transform a circle into a figure similar to a wing 
of an airplane. This figure is also looking similar to a coop profile. My idea is to transform a 
given coop profile back in order to check if the back transformation results in a “beautiful re- 
gular" figure or something like a “potato": in the latter case I would conclude that the given 
profile is not optimal in fluidic regard, it could be a “bluff body". 

Given is a circle by its radius and the coordinates of its centre. 

Applying the formulae u = x * (1 + k/(x 2 + y 2 )) and v = y * (1 - k/(x 2 + y 2 )) we obtain as trans- 
formed figure the profile of a wing in a uv-coordinate system. 

My problem is to express the variables x and y as functions of u and v. I have calculated my 
“fingers awry” but I could not get a solution. 

I was told that DERIVE should be able to solve the problem. 

It would be great if you could provide any advice or the solution of the problem. 


My (Josefs) answer was: 

First of all I’d like to present the airplane wing profile as result of the conformal mapping given in 
Mr Korner’s mail. 


#1: [,/37'COSOO - 1, V^-SINCt) + 1] 

#2: [x_ := ./37 ■COS(t) - 1, y_ := V^-SINOO + 1] 




f 25 ^ 


f 25 \ 

> 

SUBST 

x ■ 

1 + 

. y 

1 

i — i 

i i 

i — i 

l 

l 

i i 



2 2 


2 2 




^ x + y j 


^ x + y j. 

j 


#4: 


2 

2. (37. COS (t) - ./37 ■ COS (t) ■ C/37 ■ SI N (t) + 33 ) + 1 /37-SINCtJ + 32) 
2. 1 /37-C0SCtJ - 2.V37.SINOO - 39 


C2. 1 /37.SINCtJ + 2) ■ G/37 -COSO) - V37-SINOO - 7) 
2. 1 /37-C0SCtJ - 2.V37.SINOO - 39 
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I generalized the mapping and introduced a slider for k. It was interesting to observe that only for 
k= 25 a wing profile appeared. 



Next step was dealing with Mr Korner’ s problem, i.e. solving the system of nonlinear equations for x 
and y. I assumed that DERIVE - with Groebner bases implemented - should be able to do the job. 


SOLUTIONS 


#5: 


r 

r 25 i 


f 

25 ^ 

A 

c 

II 

1 1 


]_ _ 


, [x, y] 

2 2 

A V — V ■ 


2 2 




^ x + y j 



x + y ^ 

J 


# 6 : 


This did not work because the Groebner algorithm can be applied on systems of polynomials only. So 
I rewrote both equations free of fractions and then I tried again. 

I kept the system as general as possible and replaced 25 by a variable k\ 


# 7 : 


1 + 

2 2 
v x + y jj 


2 2 2 2 2 2 
(x + y ) = (u«(x + y ) = x ■ (x + y + k)) 


#S: 


v = y 


1 - 


2 2 
k. x + y JJ 


2 2 2 2 2 2 
(x + y ) = (v-(x + y ) = y-0 + y - k)) 


2 2 2 2 


2 2 2 2 


#9: S0LUTI0NS(u ■ (x + y ) = x ■ (x +y + k) a v-(x +y)=y.(x +y - k) , [x , y]) 

And I received a solution, not only one, but five including the trivial solution (pc=y = 0). But the oth- 
ers were bulky enough with two of them containing complex expressions. 

o 

2 224 224 22 2 224 2 2 4 - 

i/2 ■ ( G/CL6 ■ k + S-k-(\? - u ) + u + 2-u -v +v)-4-k+u - v ). 1 /( 1 /(16 , k + S-k-(v - u ) + u + 2-u -v + v )- 

S-u-v 


2 2 224224 22 2 2 2 4 - 

^2 , (2 , -/2 , li ■ v — C/Q-6 ■ k + S-k-(v - u ) + u + 2-u -v +v)-4-k+u - v )■ 7(70-6 ■ k + S-k-(v - u ) + u +2- 


S-u-v 


2 224224 22 2 224- 

u i/2 0/(16 ■ k + S-k»(v - u ) + u + 2-u -v +v)+4-k-u + v ) -7(7(16^ + S-k-(v - u ) + u + 2-u- 


2 S-u-v 


2 224224 22 2 224- 

u 7^ ■ i- ■ (7(16 ■ k + S-k-(v - u ) + u + 2-u-v + v)+4-k-u + v )■ 7(7(16' k + S-k-(v - u ) + u + 2-u- 


2 


S-u-v 
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I sent my results to Mr Korner and he - he does not use DERIVE - transferred the expressions to 
MS-Excel - and then had some satisfying results. He wrote (the most important parts of his letter will 
follow in a translation): 


Sehr geehrter Herr Bohm, 

nun habe ich endlich Zeit gefunden um Ihnen nochmals Dank zu sagen fur die iibermittelten 
Formeln als Ergebnis der DERIVE-Rechnungen. 

Zunachst zu Ihrer Frage woher ich Ihre Adresse bekommen habe: 

Nachdem ich vor mehreren Jahren von einem Geschaftsfreund in Siidafrika auf DERIVE 
hingewiesen wurde (er hatte eine hohere Position in einer Raffinerie und war nebenbei „Frei- 
zeitmathematiker") ist mir die Idee zur Nutzung dieses Programms fur mein anstehendes 
Problem in den Sinn gekommen. Im Internet habe ich dann bei der ..DERIVE USER GROUP" 
Ihren Namen gefunden. 

Im Internet habe ich dann auch (irgendwo) Ihren Artikel uber die moderne Unterrichtung des 
Faches Mathematik gelesen und finde ihn echt zukunftsweisend (leider finde ich ihn nicht 
mehr.... Bitte um Hilfe). 

Zu meiner Zeit und wahrscheinlich auch noch zu Ihrer Zeit wurde Mathematik sehr trocken 
gelehrt und dadurch hatte dieses Fach auch wenig Freude fur die Schuler. Ich erinnere mich 
noch an meine Oberstufe wo Kurvendiskussion besprochen, doch niemals ein Beispiel aus 
dem taglichen Leben angefiigt wurde. Dafiir wurden wir mit spharischer Trigonometrie 
(Himmelskunde) gequalt. 

Beim Studium (allgem. Maschinenbau an derTH Munchen) waren die Mathe-Vorlesungen 
auch sehr schwer und oft undurchsichtig. Ein ganzes Semester beschaftigte sich mit Funk- 
tionentheorie, konformen Abbildungen (auch 3-dimensional) und auch mit geschlossenen 
und offenen Torusformen. Im Nachhinein halte ich das fur ein Maschinenbaustudium fur vol- 
lig ubertrieben! 

Doch nun zu unserem Rechenbeispiel: 

Ich habe mir die Beschreibung des DERIVE5 - Programms von Frau Maria Koth, Universitat 
Wien, angesehen und festgestellt welche traumhaften Moglichkeiten dieses Programm bie- 
tet. 


Aus Ihren Formeln habe ich die 3. Losung und die 4. Losung programmiert und ganz eigen- 
artigerweise festgestellt, dass teilweise die Ergebnisse der 3. Losung richtig sind und teilwei- 
se die der 4. Losung. 

Damit Sie eventuell mein Beispiel mit den Zeichenmethoden von DERIVE nachempfinden 
konnen mdchte ich folgende Daten nennen: 

Der Kreis hat den Mittelpunkt M(-1;+1) und einen Radius von Wurzel aus 37. Damit kann 
man schrittweise (ich habe immer 5° genommen) die Umfangskoordinaten berechnen. 

Mit den Ihnen bekannten Formeln u = . . . und v = . . . ergibt sich dann eine Form, die einem 
T ragflugel ahnelt. 
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Die berechneten u- und v-Werte habe ich nun verwendet, um mit den langen Gleichungen 
der 3. und 4. Losung eine „Ruck“-Transformation in Richtung Kreis zu erreichen. 

Eigenartigerweise geht das dann nicht ganz ohne Probleme. Ich habe namlich festgestellt, 
dass es am Vorzeichen des jeweiligen y-Werts des Ausgangskreises liegt. Im Bereich positi- 
ver y-Werte stimmt das Ergebnis der 4. Losung, sonst das Ergebnis der 3. Losung. Beim 
Vorzeichenwechsel von y machen die Ergebnisse der 3. und der 4. Losung einen Sprung! 

Da ich aber eine bestimmte Profilform vorqeben mochte um dann mit der Rucktransformation 
zu sehen welche Form herauskommt, kann ich ja nicht vorausahnen, ob y einen Vorzei- 
chenwechsel macht Oder nicht! 

Ich hoffe doch und glaube auch, dass ich bei der Umsetzung dieser komplizierten Formeln 
nach EXCEL keinen Fehler gemacht habe!?! 

Bitte probieren Sie es vielleicht einmal mit DERIVE. 

Es gibt ja mehrere Versionen von DERIVE. Ich habe auch gehort, dass die Version DERI- 
VES besser sein soil als DERIVE6. 

Mit welcher Fassung arbeiten Sie? 

Sollten Sie mit den Derive-Plot Funktionen auch die eigenartigen Ergebnisse wie ich mit 
meinen EXCEL-Methoden erhalten, dann bitte ich Sie doch mit Mathematica nochmals einen 
Versuch zu starten. Vielleicht kommt dann irgendeine andere Losung heraus. 

Als Anhang schicke ich Ihnen meine EXCEL-Ergebnisse. Darin konnen Sie andern oderer- 
weitern wie Sie mochten. 

Fur heute viele Griifte und nochmals „Danke“. 

Klaus Korner 

Ein Teil der Excel-Datei und die zugehorigen Grafiken folgen der englischen Uber- 
setzung der wichtigen Teile des Briefes. 

. . . finally I found time to thank once more for the sent formulae which resulted from the DE- 
RIVE calculations. 

Using your expressions I took the real solutions (#3 and #4) and found out that partially re- 
sults based on #3 and partially based on #4 are correct. 

If you would like to reproduce my example with DERIVE I’ll give the following data: 

The circle is given by its centre M(-1;+1) and its radius v^7. I calculated the points of the cir- 
cumference in 5°-steps. Then using the formulae u = ... and v - ... one obtains the figure 
which is similar to an airplane wing. 

These u- and v-values were used applying the extended expressions of solutions #3 and #4 
for achieving a „back transformation “ to get the circle. 

Curiously this does not work without problems. I found out that the problems are connected 
with the sign of the respective y-values of the initial circle. For y>0 one has to take the 4 th 
solution otherwise the 3 rd one. When y changes its sign then the results of #3 and #4 show a 
jump. 
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As I want to start with a certain given wing profile in order to investigate the form resulting 
from the back transformation I cannot guess in advance whether y changes its sign or not. 

I hope - and I am pretty sure - that I didn’t make a mistake by transferring the bulky expres- 
sions to EXCEL!? 

Please try with DERIVE. There are several versions of DERIVE. I’ve heard that DERIVE 5 
should be better than version 6. Which version do you use? 


I attach my Excel-results. You are free to change or to extend the file as you like. Many 
thanks again. 

Klaus Korner 


Parts of Klaus Korner’ s Excel worksheet 


Calculations of conformal mapping and back transformation 











Centre of the circle 

V “ 

-1 

= Entry 









Vm “ 

1 

= Entry 

















Squared radius of the circle 

r 2 = 

37 

= Entry 


















Factor 


k= 

25 

= Entry 























Backtransformation 

Backtransformation 


Combination of 4th 

Deg 

Pts on Circumference 

Conformal Mapping 

4th DERIVE Solution 

3rd DERIVE Solution 


and 3rd solution 


X 

V 

u 

V 

x-new 

y-new 

x-new 

y-new 


manually combined 













0 

5,083 

1,000 

9,818 

0,068 

5,083 

1,000 

4,735 

-0,932 


5,083 

1,000 

5 

5,060 

1,530 

9,587 

0,161 

5,060 

1,530 

4,527 

-1,369 


5,060 

1,530 

10 

4,990 

2,056 

9,273 

0,292 

4,990 

2,056 

4,283 

-1,765 


4,990 

2,056 

15 

4,875 

2,574 

8,885 

0,457 

4,875 

2,574 

4,010 

-2,117 


4,875 

2,574 

20 

4,716 

3,080 

8,432 

0,653 

4,716 

3,080 

3,716 

-2,427 


4,716 

3,080 



D 


-N-L#83 


Klaus Korner: Airplane and Conformal Mapping 


P 31 







Figure 3 : Back transformation using DERIVE's 3rd solution 



Now the lower circle is the correct onel 




For 0° to 1 85° the results of the 4th solution are fitting excactly to the given circle 
The graph shows the combination of both series of data. 



One could be happy, but why do not fit the results for 355° and 360°? 

Values for 355° would be fine if we would take the 4th solution! 

1 1 1 1 1 1 — 
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I define four functions for the nontrivial solutions (sol 1 and sol 2 are the complex ones, sol 3 and 
sol 4 the real ones). 


sol3(u, v, k) := 


2 224224 22 2 2 2 4 - 

72 •((,/( 16 -k + 8-k-(v - u ) + u + 2-u -v +v)-4*k+u -v )*7(7(16-k + 8-k-(v - u ) + u + 2- 


8-u -v 


2 2 4 2 2 2 

u -v + v ) + 4-k - u + v ) + 2-72 -u -v) 


2 224 224 22 

72-7(7(16-k + 8-k-(v - u ) + u + 2-u -v +v)+4-k-u +v) v 


4 2 


sol4(u p v, k) := 

2 2 224224 22 2 2- 

72 -(2 -72 -u -v - (7(16-k + 8-k-(v -u)+u + 2-u -v +v)-4-k+u -v )-7(7(16-k + 8-k-(v 

8 • u -v 

2 4 2 2 4 2 2 

u ) + u + 2-u -v + v ) + 4-k - u + v )) v 

2 

2 224 224 22 

72*7(7(16*k +8-k-(v - u ) + u + 2-u -v +v)+4-k-u +v) 

4 

I try the functions with u and v being the parameter form of the wing - and k = 25 - and hope to come 
back to the circle: 

2 

2 ■ (37 * COS(t ) - h/ 37 * COS(t) ■ (737 * SIN(t) + 33) + 7^7 -SIN(t) + 32) 

u_ := 

2 - h/ 37 * COS(t) - 2 -h/ 37 ■ SIN(t) - 39 

2 ■ (h/ 37 ■ SIN(t) + l).(/37-COS(t) - 737.SIN(t) - 7) 

v_ := 

2 - h/ 37 ■ COS(t) - 2*737. SIN(t) - 39 

sol3(u_ p v_, 25) 

sol4(u_ p v_ p 25) 
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This is interesting and surprising as well. I don’t receive the expected circle but I receive two circles, 
one of them being the expected one. At the first view it seems to be that I am not facing Klaus 
Korner ’s sign-problem. 


But see what is happening when plotting the circles in steps: 






There are strange arcs belonging to the other circle which I don’t know how to explain. (I must admit 
that I am not an expert in conformal mappings.) 

Now I’ll deal with the second circle. It’s centre seems to lay on y = x and it passes (5,0) and (0,5). This 
is not enough to find its parameters. I didn’t know better than to use sliders for estimating the coordi- 
nates of the centre and the radius. 


Centre: (*0.73, *0.73) 

Radius «4.36 

What is the result when applying the con- 
formal mapping from above to this (sec- 
ond) circle? 



Let’s see: 
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r 


f 25 1 

SUBST 


X - 

1 + 




2 2 


K 


s. x + y > 


25 ^ll 


2 2 

x + y J. 


[x, y], [4. 36-COS(t) + 0.73, 4.36-SIN(t) + 


0.73] 


4 4 5 

0 . 01* (436 • COS(t) + 73). (3. 1828-10 -COS(t) + 3.1828*10 -SIN(t) + 2.25377-10 ) 


4 4 5 

3.1828-10 -COS(t) + 3.1828-10 -SIN(t) + 1.00377-10 

4 4 4 

0 . 01- (436 - SIN(t) + 73)- (3. 1828-10 -COS(t) + 3.1828-10 -SIN(t) - 2.4623-10 ) 


4 4 5 

3.1828-10 -COS(t) + 3.1828-10 -SIN(t) + 1.00377-10 



Obviously we receive the same profile as result of the mapping applied to two different circles (which 
are surely related in any form, but which one? Who knows?) 

(It is also possible to generate trochoids in two different ways.) 

I am much more surprised by the two little pieces of the circles. Do you have any explanation for this? 
Who is an expert? 

Josef 


Just now, finalizing this DNL, I performed the calculations with MATHEMATICA and received the 
same results - some in another form. I found that the DERIVE results are much more compact. You 
can find parts of the MATHEMATICA treatment after Mr Korner’ s next mail. 


And I had one more question which I tried to answer for myself - and for you, too, of course. 
But see first Klaus Korner’ s response: 
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Hello Mr Bohm, 

I would like to inform you about an interesting link: 

http://www.excelformeln.de/uberuns.html 

You can find numerous meaningful (and also less meaningful) EXCEL-applications (free of 
charge). Mr Klaus Kuhnlein is the main organizer of this publication. 

Finally I’d like to thank for your efforts. Your results agree completely with my Excel- 
solutions. I will try to “back transform” the profile of a real turbine blade in order to find out 
whether it results as a more or less regular figure. 

Best regards 
Klaus Korner 

These are the MATHEMATICA - calculations and plots. 


In[1 0]:= Parametric Plot [ sol3 [ a , b, 25], {t , 0, 2* Pi}, Aspect Ratio -> Automat ic ] 



0ut[1 0]= - Graphics - 
ln[1 1]:= sol4 [u , v , k ] : = 


8 u v 2 


2 -,[2 -u 2 + v 2 + 4 k + Vl 6 * v 2 + (u 2 + v 2 - 4 Jt) 2 v 3 + 4 u 2 v 2 + ^/~2 |-u 


2 w 2 'J -u 2 + v 2 + 4 k + *\j 16 Jt v 2 + (u 2 + v 2 - 4 Jt) 2 v - 8 '^2 k 'J - u 2 + v 2 + 4 



r 


1 

16 

8 v + 4 V~2" A 

j 

-u 2 + v 2 + 4 k + y 16 k v 2 + (u 2 + v 2 - 4 Jt) 2 


j 


ln[1 2]:= Parametric Plot [ so!4 [ a , b, 25], {t, 0, 2*Pi}, Aspect Ratio -> Automat ic ] 
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At first I performed the MATHEMATICA - plots using u_ and v_ from page 32. Then I asked 
MATHEMATICA to do the substitution and received quite other expressions (c and d in In[16]). 

In[1 5]:= {x*(l + 25/ (x A 2 + y A 2)), y+(l-25 / (x A 2 + y A 2))} /. {x-> Sqrt [37 ] * Cos [t ] - 1 , y -> Sqrl 


0ut[1 


5]= {(737 cos© - 1) 


1 + - 


25 


, (737 sm©+l) 


25 


(yfl 7 cos© - l) 2 + (737 sin© + l) 2 


(737 cos© - l) 2 + (737 sin© + l) 2 J 
ln[16]:= c = [ V 37 cos (t) - l] 

d = (V3 1 sin ( t) + l] 

But the resulting plot is pretty the same. 

It might be a good exercise for students proofing the identity of c and u_ (d and v_). For us it is easier, 
let’s ask DERIVE! 



c := (^37 * COS (t) - 1). 

Precision := Exact 
Notation := Rational 


25 


1 + 


2 2 
G/37-COS(t) - 1) + (737 ■ SIN(t) + 1) ) 


u_ - c = 0 

Solutions #1 and #2 are complex, back transformation applying them gives also complex expressions 
and we cannot obtain any plot. 

I wondered what the plot of the real part (and the imaginary part, too, of course) would look like and 
I was presented a big surprise: 


RE(sol l(u_ p v_ p 25)) 
IM(sol l(u_, v_ , 25)) 



The real part was again the basic wing profile! Solution #2 gave also the wing, the imaginary part is 
mirrored wrt to the origin: 
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From the eDUG-Forum: Tensor Algebra and DERIVE DOS 

Thomas Fowler: 

I posted two new versions of the Tensor Algebra file, updating the fine work done by Hans Dudler. 

The version for Derive 5 & 6 is Tensor AlgebraExtd.mth. The version for Derive 3 (last DOS release) 
is TENSORNE.MTH. I also posted an updated and expanded explanation of the operations, Ten- 
sorAlgebra.doc. The extensions are primarily for functions in use in the study of general relativity, and 
include the Riemann Curvature function K. 

Tom 

(See pages 39-41) 


Rick Nungester 

I have Derive for DOS v3.04 working on Windows XP, but it always starts full-screen instead of in a 
small terminal-like text window. Right-clicking on the shortcut and seeing its Properties shows tabs 
General, Program, Font, Memory, Screen, Misc, Compatibility, and Summary. Tve played with vari- 
ous settings, but it always starts full-screen. 

Can Derive for DOS be made to open in a small window instead of full-screen? 

I also have Derive for Windows v5.06, but use Derive for DOS v3.04 on my circa- 1994 HP 200LX 
palmtop and want the same functionality on Windows. 

Rick 


Thomas Fowler 

Until TI decides to issue a version of DERIVE for Android devices (which would be a cash cow for 
them), you can still get most of the functionality if you have your old DOS release. The last DOS re- 
lease was DERIVE 3. 1 have it running on an Android tablet under aDOSBox. It works well, but you 
really need an external keyboard. You can speed it up by going into the menu option (from the FI 
key), selecting DOSBox, and increasing the CPU cycle variable. I am running mine at 87350, about 
20x faster than the default, which naturally emulates and 8088. 

Tom 


Aleksey Tetyorko 

IMHO The last DOS release was DERIVE 4.1 1 with the extended syntax. 

Sorry, one question: can you switch to a windowed mode by pressing Alt-Enter? I cannot! My DOS 
Derives are 3.1 1 and 4.1 1. But when Tve changed the Derive Display Mode to Text (Option, Display, 
Mode:Text), I made my Derive into a window by Alt-Enter. Do your Derive and XP behave in this 
manner? 


Aleksey 
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Rick Nungester 

Thank you Aleksey, that is just what I needed. Now I can start Derive for DOS 3.04 from a Desktop 
icon directly into a window without using Alt-Enter at all. Here is what I did: 

In Derive, make Text mode the start-up setting: 


* Start Derive (no other changes before the next steps) 

* Option, Display, Mode: Text 

* Transfer, Save, State DERIVE. INI (Y, overwrite) 

* Quit 


In Windows XP, make "window” the start-up choice: 

* Desktop icon, right-click, Properties, Screen tab, Usage:Window, OK 

Now whenever I double-click the Desktop icon, Derive for DOS 3.04 starts in a text window. Plotting 
still works, drawing lines as a series of”.” text characters instead of using graphics. But then switching 
to Option, Display, Mode:Graphics switches to full-screen and good-looking graphics. If desired, then 
switching back to text mode and doing Alt-Enter will go from nice full-screen graphics back to a 
smaller text-only window. Great! 

Rick 


Aleksey Tetyorko 

For your service. You may try using DosBox ( DOS emulator for XP), and you can do it without 
switching to Derive text mode. Tve test this method some hours before now. My DosBox version is 
0.73, maybe not the last. I cannot recommend it because of slow operation. 

Aleksey 


When informing about tensors in E. Weisstein’s Concise Encyclopedia of Mathematics I came across 
this pretty quartic surface, a TANGLECUBE, Josef 



x 4 - 5x 2 + / - 5y 2 + z 4 - 5z 2 + 1 1 . 8 (1 5) = 0 
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**************************************************** 

* The Derive Tensor Algebra and Analysis Package * 

* Documentation for the utility file Tensor.mth * 

**************************************************** 

Written by Hans A. Dudler (619 420-1787) 25 November 1994 
Updated for Derive 5 1 0 November 2001 

New functions added, documentation expanded by Tom Fowler, September 201 1 

This file describes the DERIVE utility file Tensor.mth which provides definitions and functions 
for tensor algebra and analysis. Tensor.mth should normally be loaded as a utility file using 
the File > Load > Utility File command. Note that Tensor.mth activates the case sensitive 
and word input modes. As an introduction to the concepts described in this document the 
demonstration files Tensorl.dmo and Tensor2.dmo are recommended. They should be loa- 
ded using the File > Load > Demo File command. 

1. Tensor Representation in DERIVE. 

From the algebraic point-of-view tensors are an extension of the vector/matrix concept to 
higher order arrays. Thus a tensor of rank (order) 1 is a vector. A rank 2 tensor is a matrix, 
which DERIVE treats as a vector of vectors. A rank 3 tensor is a vector of matrices, a rank 4 
tensor is a matrix of matrices, and so on. The obvious way to handle tensors in DERIVE is 
therefore to extend the vector concept to higher ranks (orders). 

The equivalent of the ELEMENT function to extract an element from a tensor is the function 
(see function descriptions below) 

EL_T(A,iv) , 

where A is the name of the tensor and iv is the index 'vector'(quotes used to distinguish a 
DERIVE vector from a rank 1 tensor). As an example consider the tensor of rank (order) 3 
and dimension 3 

a b c 

A= d e f 
[g h i 

Then 

EL_T(A,[2,3,2]) returns q 
EL_T(A,[3,1,3]) returns u 

As is clear from the example, the first index (element 1 of the index vector) selects the ma- 
trix, the second index selects the row-vector in the selected matrix and the third index finally 
selects the element in the row vector. Note that under DERIVE the indices must range from 1 
to n (the space dimension), zero (often used in relativity) is not allowed. From the above ex- 
ample we deduce that 

E L_T ( A, [i , j , k] ) <=> ELEMENT(ELEMENT(ELEMENT(A,i),j),k) 
in this case, and the extension to higher ranks is obvious. Note that 
EL_T(A,[i]) <=> ELEMENT(A,i) 



First page of the documentation, one example follows: 
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In this case we have chosen the last of the three possibilities for t, and assumed that a>0, 
giving t = ATAN(s/a), and converting T to s variable by substituting this and invoking Simplify, 
Basic. This allows us to get rid of the SIGN function, so we have 



2 2 

' s y 

#110: T := 

V(a + s ), ATAN 




s, 3 ) . 


Use these values for xi’s in terms of s to substitute into G : 


#111: G_ := 


0 0 



1 


2 

+ s ) 


1 



2 2 

- ./(a + s ) 0 


Invoking CURVATURE on this will yield 0, which indicates that the original curve was a geo- 
desic. 

#112: CURVE_T(T, G_, s) = [0, 0] 


(There are problems with the ‘T” in the DERIVE-file because T is reserved for the Gamma-function. I 
took G_ instead, Josef) 


About Tensors (from CRC Concise Encyclopedia of Mathematics): 

An nth-rank tensor is a mathematical object in m-dimensional space which has n indices and 
m n components and obeys certain transformation rules. Each index of a tensor ranges over 
the number of dimensions of space. If the components of any tensor of any rank vanish in 
one particular coordinate system, they vanish in all coordinate systems. 

Tensors of rank 0 are calles scalars, tensors of rank one are vectors and tensors of rank two 
are matrices, ... 


Two recommended links to introductory papers on Tensors 

An Introduction to Tensors for Students of Physics and Engineering (NASA): 
http://www.qrc.nasa.qov/WWW/K- 

12/Numbers/Math/documents/Tensors TM200221 1716.pdf 


A paper in GERMAN (Technische Universitat Graz): 

http://www.math.tuqraz.at/~qanster/lv vektoranalvsis ss 10/17 tensor definition.pdf 
(GERMAN) 
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Von: R.N. Folsom [rnfolsom@redshift.com] 

Gesendet: Donnerstag, 16. Dezember 2010 09:24 
An: "Bohm, Josef" 

Betreff: A Constrained Inequality's Truth 

Josef : 

First, thank you very much for your message of 13 December. I 
very much appreciate your work supporting the Derive User Group, and 
me . 

I have attached a zipped file, "GrowthRates Summed vs Multi- 
plied, N3 - Derive, " that contains a dfw file of the same name, 
which describes a problem I am having. I hope you can pass the file 
on to someone who might know the answer, or that you decide it is 
worth including in the next newsletter bulletin. 

But I would be happy to revise this document if you think I 
should do so. 

Sending the Comparative Statics Technique information is defi- 
nitely on my short term todo list. 

Thanks again for your help. 


As you will see below, I am trying to formally prove the circumstances under which a particu- 
lar inequality is true — and Derive's Simplify or Solve command actually says that it is true. 
(For now, all I need is sufficient conditions.) 

But I have not been able to figure out the Derive6.1 command or commands that will do that. 

My sources are Bernhard Kutzler and Vlasta Kolol-Voljc's "Introduction to Derive 6, Ad- 
vanced Mathematics for Your PC" (hereafter Derive6); and Albert and Joan Rich, Theresa 
Shelby, and David Stoutemeyer's "Derive User Manual Version 6" (hereafter Derive3). 

Derive6's equation solver (page 57), when used to determine whether a relationship is True 
or False, doesn't apply to what I am trying to do because it doesn't allow conditions to be 
imposed on the equation's (or, in my case, the inequality's) variables. 

Derive6's discussion of the IF tool (on pages 166 and 227) didn't help because I didn't un- 
derstand page 166 (I don't know any Analytic Geometry), and although I think that I did un- 
derstand page 227 I couldn't figure out how to apply it to my problem. 

Derive3's discussion of the IF tool (pages 309-310), made sense, but again I couldn't figure 
out how to apply it to my problem. 

My understanding is that the IF tool has four fields: a test condition, a "then" clause of re- 
sults if the test condition is satisfied, an "else" clause of results if the test condition is NOT 
satisfied, and an "unknown" clause. 

Derive3's opening IF example (page 309) is a payroll problem: IF(h <= 40, 1 0h, 400 + 15(h- 
40)) where h is hours worked and lOh and 15(h-40) are worker payment amounts. Unfortu- 
nately, the problem does not include the unknown clause, which if not included but is needed 
turns out to be the entire IF(...) tool. An example that does include an unknown clause would 
have been helpful. 
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In statement #6, multiplying out the (rjn + 1 ) factors would eliminate the -1 - (r|3 + r|2 + rjl ) 
expression. Doing that (by using Derive's Simplify-Factor command) gives statement #7. 


-1 < r .(r -Cr + 1) + r ) + r ■ r 
#7: 321 1 12 


or 


-1 < r ■ r ■ r + r • r + r • r + r ■ r 
#8: 128 28 18 12 

Note: I have done the same calculations above using N = 5. 

Given the sufficient condition that all rjn > 0, statements #7 and #8 are obviously true be- 
cause the left side is negative, and every term on the right side is positive. 

And since statements #7 and #8 are merely a rearrangement of #5, which is an example of 
statements #4, #3, and #2, statement #7 demonstrates that if all r|n > 0, the "usual situation" 
inequality holds, so that summing N periodic rates of change rjn (instead of multiplying cor- 
responding 1 + rjn numbers) understates the true total rate of change over N periods. 

HOWEVER, I need to learn how to persuade Derive6.1 to do a formal proof (results either 
true, false, or unknown) of the preceding paragraph's results, because my next step is to 
experiment with less restrictive conditions, such as all r|n >= -1. 


TO DO: Figure out a formal proof — for any statement #5 through #8 — that when #10 (be- 
low) holds, comes up with "true" rather than the wrong answer of "false." 

For example, Simplify-Basic of either IF(#10, #7) or IF(#7, #10) returns the same (written out) 
IF statement, presumably meaning that Simplify could not determine whether the inequality 
was true or false, despite the "all rjn >= 0" variable conditions. 

But Solve-Expression of either IF(#10, #7) or IF(#7, #10) concludes that the constrained ine- 
quality is false, when I know that it is true (from the analysis that follows #8). 

I clearly am not using IF properly. It may not be intended for what I am trying to do. In that 
case, I need to know what else to do. 

VECTOR (r >0, n, 1, 3, 1) 

#9: n 

|"r > 0 a r > 0 a r >01 

#10: L 3 2 1 J 

Roger Folsom 

Thanks for welcoming me to the Derive Users Group. 

Actually, it turns out (prior to the Derive 6.1 CD) that al- 
though my latest MS-DOS diskettes were for Derive 4.0.7, I have a 
Derive 3 manual (in addition to the Derive 6 manual by Bernhard 
Kutzler and Vlasta Kokol-Vol j c) , so apparently my use of Derive goes 
back at least to Derive 3.x. The Derive 3 manual is quite useful in 
providing more detail about various Derive features and capabili- 
ties . 


> Applications in Economics would be highly appreciated! ! < 
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You may already have some. I have, on paper, a message posted 
by "R.N. Folsom” on the Derive BBS, 26 April 1993, about "Folsom’s 
Comparative Statics Technique." I was responding to a message 
#2896, 19 April 1993, on "Derive/Mathware" (apparently a different 
BBS) by "Harald Lang to Roger Folsom," with the title "Harald Lang's 
Use of Taylor's Series to Specify the Chain Rule." 

And that was preceded, again on "Derive/Mathware," by message #2893 
by "Hadud to Roger Folsom about Response To Questions," which appar- 
ently was a response to a set of 7 questions posted by me. Unfortu- 
nately, I don't have the seven questions, but Hadud' s answers give 
pretty good clues about what my questions were. However, they ap- 
parently were about how Derive worked, and did not have anything to 
do specifically with economics. 

But "Comparative Statics" is something that economists do fre- 
quently. If you have a system of functions i.e., a model 

describing a set of markets (as in microeconomics) or describing in 
broad terms the entire economy (macroeconomics) , the variables fall 
into two categories: endogenous (whose values are determined by the 

model's functions) and exogenous (whose values are determined out- 
side of the model, e.g. Russia's awful wheat-destroying fires this 
past summer) . And you want the model to tell you what the exogenous 
event will do to markets (not merely the market for wheat, but also 
to markets for grain substitutes such as corn, other things that 
grow such as cotton, markets for beef since cattle often are fed 
corn, etc.). in short, you want to know (more accurately, to esti- 
mate) the effects of an exogenous variable (or possibly multiple ex- 
ogenous variables) on endogenous variables such as the prices and 
quantities of corn, cotton, beef, etc. 

Comparative statics usually (maybe always) involves the "chain 
rule": derivatives (or finite difference quotients) of y = Y (x) and 

x = X(a) where x and y are endogenous and a is exogenous, and where 
I have used capital letters instead of "f" to indicate a function 
(because I'll probably want to use f for something else, such as 
fire severity!) . In short, you want to know not merely dy/dx and 
dx/da but also dy/da = [dY (x) /dx] . [dX (a) /da] , which in the U.S. and 
maybe elsewhere is commonly known as the "chain rule." 

But for a multi-function system, to calculate a full set of to- 
tal derivatives of all of the endogenous variables with respect to 
each of the exogenous variables, Derive wouldn't use the chain rule. 
Harald Lang came up with one solution, and I used his notation to 
come up with another solution (namely, not substituting any function 
into any other function before taking derivatives and calculating 
Jacobian matrices for the system) . 

If you do not already have the messages I have described above, 
let me know and I will send them. I'd do so right now except that 
they all are paper documents and the ink is not as dark as I wish it 
were, so I would need some time to fuddle around with my scanner and 
optical character recognition software to send them electronically. 

I don't have time to do that right now because I am facing deadlines 
for two entirely different projects, but I think I could get to it 
in January. 

My apologies for taking your time in trying to "explain" the 
terms comparative statics, and endogenous and exogenous. I did so 
because I don't know if those terms are used by non-economist mathe- 
maticians, or in countries outside the U.S. and Canada. 
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I briefly visited your Didactics of Computer Algebra 
http : //www . acdca . ac . at/ , and tentatively concluded that it focuses 
on teaching computer algebra. But the site's index frequently uses 
the CAS acronym, where presumably C is for Computer and A is for Al- 
gebra, but I'm wondering what the S is for. Studies? Or a non- 
English word? At your convenience, please satisfy my curiosity! 


A couple of suggestions: After downloading all of the zipped 

MTH files from the DUG website, one at a time, I then saw and fol- 
lowed the recommendation to look at the "back" of the PDF membership 
file, which contained the statement "Disk of the year containing all 
.MTH files is included." It's been a long day (almost 2:00 a.m. 
California time), so I likely am operating in stupid mode, but: In 

what is that "Disk" included? 


With or without a disk, it might be useful if the DUG website had 
three large zip files that included each of the three groups of num- 
bered MTH zip files. 

Similarly, since the newsletters/bulletins apparently include arti- 
cles that explain the use of those zipped mth files, it might be 
useful if the DUG website had zipped packages of dnl*.pdf files, so 
each pdf file wouldn't have to be downloaded individually. 

I'll close with a very minor request: Mathematicians, and also 

Derive, like to write polynomial functions "downhill, " meaning with 
the largest exponent term first, followed by successively smaller 
exponent terms, and then a constant. I do understand that that se- 
quence has mathematical advantages, for example by facilitating di- 
viding one polynomial by another. But economists often like to 
write polynomials "uphill" (constant term first, followed by terms 
with successively larger exponents), because we often don't know how 
big the polynomial needs to be until we use data to estimate its co- 
efficients. And it is more convenient to add or delete a high-order 
term if it is at the end of the polynomial rather than at its begin- 
ning, because if you add or delete a term at the beginning of the 
polynomial the coefficients of the previous version (if identified 
by subscripts) have to be renumbered. 

I wouldn't want to change the way Derive works. But if you know of 
a Derive operation (probably created by some user, since I'm pretty 
sure it's not in Derive itself) that rearranges a downhill polyno- 
mial (and any polynomial subcomponents) and rewrites it (and its 
subcomponents) to go uphill rather than downhill, please let me 
know. That would be a lot faster than doing reordering the polyno- 
mial manually, particularly if it is in a factored form. 

Thanks again for welcoming me into the Derive Users Group. 
Cordially, Roger Folsom 


BTW, it is not easy to find a meaningful translation of “Comparative Statics” into German, Josef 
Interesting readings on this issue can be found among others at: 

http://en.wikipedia.org/wiki/Comparative statics 
http://www.applet-magic.com/compstat.htm 
http://www.economics.utoronto.ca/osborne/MathTutorial/DIFF.HTM 
http://www.duke.edu/~dgraham/handouts/ComparativeStatics.pdf 
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Some days ago I was phoned by Walter Klinger, a DERIVIAN from the first days, about an interesting 
behavior of DERIVE. It must be said that Walter has used several technologies in math teaching 
(Lower and Upper Secondary) during the last years and finally came back now to his “first love”, 
DERIVE. He explained the problem on the mobile phone and then sent a short DERIVE file to explain 
his problem more detailed: 


An Equation of Equations 

Walter Klinger, Stockerau, Austria 


Textbook 5th form (15 years), Example 899: 


#1: [g := 5-x - 12-y = -3, h := 5-x - 12-y =5, n := [5, -12]] 

#2: g = h 

??? What do we expect when solving this equation? 
Input produced by students ! 

In fact there are two equations with tow unknowns! 

#3: SOLVE (g = h, [x , y]) 


#4: 5-x - 12-y * 5 a 5-x - 12-y i -3 

??? Strange answer? How to interpret this answer? 
Lets plot it? 

What is an equation of two equations? 

Plot of expression #2 
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My first answer via phone was, that I cold imagine that DERIVE is interpreting h and g as statements 
and then treats expression h = g as a Boolean expression. 

But I was not really sure and asked myself for an explanation for #4. 

What to do? Simply ask DERIVE applying its really unique Stepwise simplification'. 

#7: g = h 

#8: 5 -x - 12-y = -3 = (5-x - 12-y = 5) 

If p and q are of logical type, 

p = q -» p « q 

#9: 5-x - 12-y = -3 « 5-x - 12-y = 5 

p «• q -> (p ■+ q) a (q ■+ p) 

#10: CCS -x - 12-y 5-x - 12-y = 5) a (5-x - 12-y = 5 -» 5-x - 12-y)) = (-3 « 5-x - 12-y = 5) 

p^q->-pvq 

#11: ((- 5-x - 12-y v 5-x - 12-y = 5) a (5-x - 12-y = 5 -» 5-x - 12-y)) = (-3 « 5-x - 12-y = 5) 

The first steps confirmed my conjecture about treating g = h as a Boolean expression. Some steps are 

following until coming to an end which is given in #16: 


If n and m are integers, 

n a m — * bitwise -logical -and-of -n-and-m 

#16: ((5-x - 12-y = 5 v 5-x - 12-y) a (5-x - 12-y # 5 v 5-x - 12-y)) = ((5-x - 12-y = 5 v 2) a 

(5-x - 12-y * 5 v -3)) 


Repeating ... 

This would be the correct form of the Boolean expression 
#17: q a h 

which gives the correct truth value for the contradiction 

#18: false 

My explanation for #4 from above 

#19: 5-x - 12. y * 5 a 5-x - 12-y ± -3 

The expression g = h or better g <-> h is true if 5-x - 12-y 1 5 AND 5-x - 12-y i- -3 
i.e. it is false for 5x- 12 = 5 OR 5x - 12y = -3. 

Hence the whole plane is shaded (true) with exeption of (only some) the points on the lines g and h. 
Why only some points? This could be a deficiency of the shading algorithm?? 

I tried with two other equations of parallel lines and found out that the “non-plotting” of the 
“forbidden” points seems in fact to depend on the resolution of the screen and connected with this on 
the scaling of the plane. 


See the next example: 
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#1: [h 3-x + 5-y = 12, g := 3-x 

#2: SOLVECh = g) 

#3: 3-x + 5-y t 12 a 3-x + 5-y t 

Many white spots within the 


+ 5-y - -5, j := 2-x + 5-y 
-3 

shaded plane 


5 ] 


















' . ' . . ' . ' . ' . 
















: . ' ' . ' ' . 


















































I change the scaling and there are no spots at all! 
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There are the spots again but not the same as above 



My next question: What will happen in case of a solution? Take equations h and j. 

#4: SOLVE (h = j) 

#5: Ox + 5-y = 12 v 2-x + 5-y t S') a (2-x + 5-y=5v 3-x + 5-y t 12) 

The (hpj) system has a unique solution, 
but ... #7 cannot be plotted 

#6: S0LVE(h = j, x) 

#7: 


' 12 - 5-y 5-(l - y) ^ 

V A 


r 12 - 5-y 5.(1 - y) ^ 

-f- A -f- 

.■ 1 jS 

Lb 2 J 

V 

j\ -f- A -p— 

Lb 2 J 


SOLVE (h = h) 
#9: true 


Any comments on this interesting behavior are welcome. I tried other Computer Algebra Systems. 
Most of them do not accept the “equation of equations”. 

I tried with changed resolution of the screen (1024 x 768 instead of 1280 x 1024) and received other 
patterns in the shaded plane. 


